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Abstract. We consider the aggregation equation Ut = V' (Vu — uVK.(u)) in a bounded 
domain f2 C R d with supplemented the Neumann boundary condition and with a non- 
negative, integrable initial datum. Here, K. = JC(u) is an integral operator. We study the 
local and global existence of solutions and we derive conditions which lead us to cither 

en 

, the stability or instability of constant solutions. 

G\ \ 

I. Introduction 

^ . We consider the initial value problem for the following non-local transport equation 

-S ■ (1.1) u t = V ■ (Vu - uVKiu)) for x e Q C M d , t > 0, 

2 ■ 

with supplemented the Neumann boundary conditions i.e 

(N ■ (i.2) — = o for x e an, t > 

j> ! on 

m 

On 
(N , 

in. (1.3) u(x,0) = tto(s). 

O ■ Here, the operator K{u) = )C(u)(x,t) depends linearly on u via the following integral 



X 



and a nonnegative initial datum 



formula 

(1.4) K(u)(x,t)= I K(x,y)u(y : t)dy 



for a certain function K = K(x, y) which we call as an aggregation kernel. 
There is large number of works considering the inviscid aggregation equation 

(1.5) u t + V- (u{VK*u)) = 
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in the whole space M. d which has been used to describe aggregation phenomena in the 
modelling of animal collective behaviour as well as in some problems in mechanics of 
continuous media, for instance, [51 EHl EE]. The unknown function u = u(x,t) > repre- 
sents either the population density of a species or, in the case of materials applications, a 
particle density. Equation (II. 5p was derived from the system of ODE called "individual 
cell-based model" [6j [22] representing behaviour of a collection of self-interacting particles 
via pairwise potential which is describe by aggregation kernel K. More precisely, equation 
(11. 5p is a continuum limit for a system of particles Xk(t) placed at the point k in time t 
and evolving by the system of differential equations: 
dX k (t) 



VK(x k (t)-Xi(t)), kez 



dt 

iez\{k} 

where K is the potential. 

Questions on the global-in-time well-posedness, finite and infinite time blowups, as- 
ymptotic behaviour of solutions to equation (ll.5p . as well as to the equation with an 
additional diffusion term, have been extensively studied by a number of authors; see e.g. 
[H El El El E2 H3l E] and reference therein. 

One introduces the diffusion term in (jl.5p to make the model more realistic and to 
describe the interesting biological (and mathematical as well) phenomenon: competition 
between aggregation and diffusion, see e.g. [U El EH EB] . 

In this work, however, our main motivation to study such models is that, in particular 
case, equation (11.11) corresponds to the parabolic-elliptic system describing chemotaxis, 
namely: 

(1.6) u t = V • (Vu-uVv), -Av + av = u, i£fl, t > 

for a positive constant a. In this system, the function u = u(x,t) represents the cell 
density and v = v(x,t) is the concentration of the chemical attractant which induces a 
drift force. Here, the function K (x, y) is the Green function of the operator — d 2 x + al on 
Q with the Neumann boundary conditions. Moreover, it is called the Bessel potential and 
it is singular at the origin if d > 2. On the other hand, in the one-dimensional case, when 
Q = [0, 1] and a = 1 this fundamental solution is given by the explicit formula i.e. 

1 pX+y _i_ p 2-x-y I x-y _i_ y-x 

(1.7) K( x ,y) = v -™+ e + % (e2 ^ 

In this work we derive some properties of solutions of aggregation equation in a bounded 
domain under no flux boundary condition (11.21) . The main goal, is to study stability of 
constant solution. In particular, we derive conditions under which constant solutions to 
problem (ll.ip - (ll.3l) are either stable or unstable. Here, let us point out that instability 
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result does not depend on dimension of the domain, and cover the case when the ag- 
gregation kernel comes from chemotaxis model (11.61) . Hence, even though solutions are 
global-in-time and bounded, a constant steady state can be unstable. This mean that 
even in one-dimensional chemotaxis we can observe the competition between aggregation 
and diffusion mentioned above. 

For the completeness of exposition we also discuss existence of solutions to (jl.ip - (11.3p . 
In order to do that, we use techniques which are rather standard and well known. In 
particular, we show that under some general condition on aggregation kernel we can al- 
ways construct local-in-time solution to (ll.ip -f Tl~3|) . However, some additional regularity 
assumption on the initial datum have to be imposed if V X K is in some sense too singular. 
Moreover, for mildly singular kernels (see Definition 12.91 for precise statement), problem 
( ll.ip -( fTT3"j) has a global-in-time solution for any nonnegative and integrable initial condi- 
tion. 

Notation. In this work, the usual norm of the Lebesgue space L P (Q) with respect to the 
spatial variable is denoted by || • || p for any p e [1, oo] and W k,p (Q) is the corresponding 
Sobolev space. The letter C corresponds to a generic constants (always independent of 
x and t) which may vary from line to line. Sometimes, we write, e.g. C = C(a, /3,7, ...) 
when we want to emphasise the dependence of C on parameters a, (3, 7, .... 

2. Main results and comments. 

2.1. Stability and instability of constant solutions. In this paper, we assume the 
following conditions on the aggregation kernel 

dK 

(2.1) _(. ;2/ ) = onffi for all y G Q, 

(2.2) V, / K(x,y)dy = 0, 

Jn 



(2.3) ||V x .if||oo, g ' = ess sup \\V x K(x, -)\\g> + ess sup \\V x K(-,y)\\ q > < 00 

for some q E [1, 00]. 

Remark 2.1. Notice that under the assumptions (12. ip . a solution of problem (ll.ll) - (ll.3l) 
conserves the integral (the "mass") i.e. 

(2.4) ||n(t)||i = / u(x,t)dx= / u (x) dx — \\u Q \\i for all t > 0. 

Jn Jn 
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Indeed, it is sufficient to integrate the equation ( II. ip with respect to x and use identities 
( II. 2p and (12. ip . Moreover, this solution remains nonnegative if the initial condition is so, 
due to the maximum principle. 

Remark 2.2. Note, that assumption (12. 2 p implies that every constant function u = M 
satisfies equation ( II. ip . In fact, the chemotaxis model ( II .6p is our main motivation to 
state this assumption. Indeed, if (U, V) is a stationary solution to (II. 6p then U is constant 
if and only if V is constant as well. It means that, if the kernel K is the Green function 
of the operator —A + al then the term VfC(U) in equation (ll.ip for U = M, has to be 
equal and so, K satisfies (12. 2p . 

The main goal of this work is to study stability of constant solution to problem ( II. 1 p — 
( 1 1 . 3 p . More precisely, we look for sufficient conditions either on the stability of constant 
solutions or their instability. Our result can be summarise in the following way 

• If the constant solution u(x,t) = M > of problem (ll.ip — (11.3P is sufficiently 
small, then it is asymptotically stable solution in the linear and nonlinear sense, 
see Proposition 12.31 and Theorem 12.41 below. 

• If the constant solution u(x,t) = M > is sufficiently large, then there is a large 
class of aggregation kernels (which include the kernel coming from chemotaxis 
system (11.61) ). such that u(x,t) = M is a linearly unstable solution of (ll.ip - (ll.3p . 

Thus, we focus on a solution to problem ( ll.ip — (11.3P in the form 



where M is an arbitrary constant and tp is a perturbation. Moreover, we assume that 
J Q (p(x, t) dx = for all t > 0, to have 



Hence, from equation (ll.ip . using assumption ( I2.2p . we obtain the following initial bound- 
ary value problem for the perturbation tp 



u(x,t) = M + <p(x,t), 




(2.5) 




(2.6) 



TT- = for x G dQ, t > 
on 

<p(x,0) = (p (x). 



(2.7) 
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We also introduce its linearized counterpart, namely, we skip the term V • (</?V/C(y?)) 
on the right hand side of (12. 5p to obtain 

(2.8) <p t = A(p-V- (MVK{<p) 

(2.9) ^ = for x G dQ, t > 

on 

(2.10) <p{x,0) = <p o (x). 



In the following, we use the linear operator Cip = —Atp + V • ^MV/C(<£>)J with the 
Neumann boundary conditions, defined via its associated bilinear form 



(2.11) J(y?,^)= / V^-V^dx-M / V/C(^)V^dx 

Jn Jn 

for all <p,i/j G W 1 ' 2 ^). 

Here, we recall that a constant M is called a linearily asymptotically stable stationary 
solution to nonlinear problem (II. II) - (II. 3p if the zero solution is an asymptotically stable 
solution of the linearized problem (I2.8l) - (l2.10p . Moreover, a constant M is called linearily 
unstable stationary solution to nonlinear problem (ll.ll) - (ll.3p if zero is an unstable solution 
to linearized problem (12T81- (I2~TUD . 

Proposition 2.3 (Linear stability of constant solutions). Assume, that the aggregation 
function K(x,y) satisfy conditions (12. ip and (12. 2p . If, moreover, the operator VIC : 
L 2 {VL) — > L 2 (f2) given by the form V/C(</?) = f n V x K(x,y)(p(y) dy is bounded and if 



(2.12) M||V/C|| L2 ^ L2 < y/\ h 

where Ai is the first non-zero eigenvalue of —A on Q under the Neumann boundary condi- 
tion then M is a linearily asymptotically stable stationary solution to problem (ll.ip - (ll.3p . 

We prove this proposition in Section |3j Here, we only emphasise that proof allow us to 
show the nonlinear stability of constant steady states. Under slightly stronger assumptions 
imposed on the kernel K. 

Theorem 2.4 (Nonlinear stability of constant solution). Let the assumptions of Propo- 
sition \2.3\ hold true. If moreover \\ VxKW^^ < 00 then there exists a positive con- 
stant n = n(V x K, M,Q) such that for every ipo G L 2 (Q) satisfying H^olh < V an d 
j n ipo(x) dx = 0, the perturbed problem (I2.5p - (l2.7p has a solution ip G C([0, oo), L 2 (Q)) 
such that J n ip(x,t) dx = for allt > 0. Moreover, we have 

\\ip(t)\\ 2 -> as t ->• oo. 



Next, we discuss instability of constant solutions. 
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Theorem 2.5 (Instability of constant solutions). Let w\ = W\(x) be the eigenf unction 
of —A on Q under the Neumann boundary condition corresponding to the first nonzero 
eigenvalue Ai, and such that || 1 1| 2 = 1- Assume that ||V/C||i2^ L 2 < 00. If moreover, the 
aggregation function K(x, y) satisfy 

(2.13) / / K(x, y)w 1 (y)w 1 (x) dx dy = A > 0, 

Jn Jn 

then for M > 1/A the constant solution M of problem (II. ip - 01.31) is linearily unstable 
stationary solution. 

Remark 2.6. Let us notice that the aggregation function K which comes from chemotaxis 
model (jl.6p satisfies the condition f)2.13p . Indeed, in this case, K(x,y) is a fundamental 
solution of the operator —A + ai in a bounded domain supplemented with the Neumann 
boundary conditions. Thus, the function 



w(x) = / K(x,y)w 1 (y) dy 
Jn 

satisfies the following equation 

(2.14) -Aw + aw = w 1 . 

After multiplying equation f!2.14p by W\ and integrating over Q and using the Neumann 
boundary condition we obtain 

Aoti dx + a wwi dx = (wi) 2 dx. 
n Jn Jn 

Obviously, by the definition of A, we have J n ww\ dx = A. Thus, after integrating by 
parts we obtain 



(2.15) - / wAw l dx = 1 - aA. 

Jn 

Finally, we use the fact that wi is the eigenfunction of —A to get 

Ai / wwi dx = 1 — aA, 
Jn 

which implies that A = > q 

Remark 2.7. Our stability results on constant steady states corresponds to the well-known 
results on the global existence versus blow-up of solutions to Keller-Segel system (11.61) . 
In particular, for general kernels (see Definition 12.91 below), where solutions of problem 
fll.ip - fll.3p are global-in-time, we can still observe the competition between the diffusion 
and the aggregation. 
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Remark 2.8. Let us mention, that the inviscid aggregation equation ( II. 5p in the whole 
space M. d can be formally considered as a gradient flow of the energy functional 

E( u ) = 77 / / K(x — y)u(x)u(y) dx dy 

with respect to the Euclidean Wasserstein distance as introduced in [2U] and generalized to 
a large class of PDEs in [8] and in [7]. We have proved that, in some sense, if this energy 
functional on the first eigenfunction of —A is positive then sufficiently large constant 
solutions of the system (1.1)-(1.3) are unstable. 

The proofs of Theorems 12.41 and 12.51 are given in Section [3j 

2.2. Existence of solutions. For the completeness of exposition we also study the exis- 
tence of solution to f ll.ip - fll.3p . First, let us introduce terminology analogous to that one 
in [T2]. 

Definition 2.9. The aggregation kernel K : Q x Q — y IR is called 

• mildly singular if || V x K\\oo,q' < oo for some q' G (d, oo]; 

• strongly singular if || V^i^Hoo^' < oo for some q' G [l,d] and \\ V^i^Hoo^' = oo for 
every q' > d. 

Notice that aggregation kernel taken from one dimensional chemotaxis model fll.6p is 
mildly singular in the sense stated above. 

We begin our study of properties of solutions to the initial value problem f ll.ip - fll.3p 
by showing the existence of solutions which depends on the quantity || Va^Hoo^' defined 
in ( 12~3|) . 

First, we show that for mildly singular kernels, solutions to the problem fll.ip - fll.3p are 
global in time. 

Theorem 2.10 (Global existence for mildly singular kernels). Assume that there exists 
q' G (d, oo] such that \\V x K\\oo :q ' < oo where ||V :E ii'|| 00ig / is defined in f)2.3p . Denote 
q = ^ [i 5 g]/ (d — 1)). Then for every initial condition uq G such that Uq{x) > 

and for every T > problem fll.ip - fll.3p has a unique mild solution in the space 

y T = C{[0,T],L 1 (Q))n{u:C([0,T],L q {Q)), sup t'^IMI, < oo} 

0<i<T 

equipped with the norm \\u\\y T = sup 0<t<r ||w||i + sup 0<t<T t 2 ^ 



Next, we show the local- in-time existence of solutions to fll.ip -f TO]) for the case of 
strongly singular kernels. 
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Theorem 2.11 (Local existence for strongly singular kernels). Assume that there exists 
q' G [1, d] such that || V x i^||oo, g ' < oo. Let q G [d/(d — 1), oo] satisfy 1/q + 1/q' = 1. Then 
for every positive u G R L q (Q) there exists T = T(||tt ||i, ||wo||g> ||V x .K'|| 00| g') > 

and a unique mild solution of problem (jTTTJ) — (JTT3J) in the space 



equipped with the norm \\u\\x T = sup 0<t<T ||u||i + sup 0<i<T \\u\\ q . 

Remark 2.12. Let us mention, that our previous, stability results imply the global-in time 
existence of solutions for strongly singular kernels provided initial data are sufficiently 
small. More results on the global-in-time solutions to fl 1 . 1 [) - ( TTT3|) in the whole space 
Q = M. d can be found e.g. in [L2] . 

Remark 2.13. Karch and Suzuki in their work [12] studied the viscous aggregation equa- 
tion, namely the equation (II. ip considered in the whole space M. d . They show that there 
are strongly singular kernels (in the sense similar to Definition 12. 9p . such that some so- 
lutions blow up in finite time. Moreover, there is a large number of works studying the 
blow-up of solution to chemotaxis model (jl.6p . see e.g. [T71 Q21 [TH] and reference therein 
as well as the review paper by Horstmann [10] for additional references. 



eigenvalue of —A on Q under the Neumann boundary condition. 
Now, we are in the position to prove the Theorem 12.31 

Proof of Proposition \2.3[ After multiplying equation ( 12. 8 p by ip and integrating over Q 
we get 



Xr = C([0, T], L\Q)) n C([0, T], L"(n)) 



3. Stability and instability of constant solutions 



In our reasoning, we use the following Poincare inequality 



(3.1) 




which is valid for all ip G W 1,2 (Q) satisfying j n tpdx = 0, where Ai is the first non-zero 




Now, using the Cauchy inequality we obtain 



(3.2) 




STABILITY OF SOLUTIONS TO AGGREGATION EQUATION 9 

Finally, we apply the Poincare inequality (13.1 j) to get the following differential inequality 

d 



±y(. it )\\j<(-\ 1 + M 2 \\vic\ 



L*-H.*)\m\2 



which, under assumption (|2.12p . directly leads us to the exponential decay of ||</?(t)||2 as 
t ->• oo. □ 

Proof of Theorem 2.4 ■ After multiplying equation ( 12. 5 p by ip and integrating over Q we 
get 

(3.3) - — = -J(<p t <p) + / ^V/C(^)V^dx 

where J is the bilinear form defined in (12. lip . In ( 13. 2 p , we have already got the inequality 

1 f M 2 f 

(3.4) -■%,¥>) < ~2 / |V^| 2 dx + — ||V/C||i 2 ^ L2 / ^ 2 dx. 

To estimate the second (nonlinear) term on the right-hand side of (I3.3p . we use the 
e-Cauchy inequality, as follows 

/ <^V/C(^)V^dx < e [ (V^) 2 dx+^- f </? 2 (V/C(y?)) 2 dx 
Jn Jn 4e J n 

(3.5) <£ I (V^dx+^HV/CMUL / ^ 2 dx 

^(V^ + ^^/^d*) 2 , 

since 

|| / V x K(-,y)ip(y) dt/Hoo < esssup \\V x K(x, OlhlMh = || V^Hoo^Hvlh- 
Jn sen 

Applying inequalities (13 ,4p and ( 13. 5 p in (13. 3p we obtain 
J if 2 dx < ( - 1 + 2e) y (V</?) 2 dx 

v^||L 2 / /• , , \ 2 



+ (M 2 ||V/C||i^ L2 ) / y 2 dx + g^^ ( I V 2 dx 



and finally using Poincare inequality ( 13. ip we get the following differential inequality 

^|M| 2 < (A, (2e - 1) + M 2 || VfC\\l^ L2 ) Ml + " V ^" L - 2 Mi 

Notice, that under assumption ( I2.12p . we can find e > small enough that the term 
Ai(2e - 1) + M 2 ||V/C|| 2 2^ L 2 ) is negative. Thus, the proof is complete because every 
nonnegative solution of the differential inequality /' < —C\f + C 2 f 2 with f(t) = \\f(t)\\l 
and with positive constants C±, C 2 decays exponentially to zero, provided /(0) is suffi- 
ciently small. □ 
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To study the instability of constant solutions, first, we consider eigenvalues of the 
operator £ defined via its bilinear form (12. lip . 

Lemma 3.1. Let the operator 

(3.6) dp = -Aip + V • (mV/C(» 

supplemented with the Neumann boundary condition be defined by the associated bilinear 
form J(ip,ip) given in (12. lip on W 1,2 (Q). Assume that V X K G L 2 {Vt x Q) satisfies (12. ip . 
Then, the number 

(3-7) A= inf 

(psw i2 (n) (Mil 

is finite and there exists (p G W /1,2 (fi) swc/i i/iai 

A ~ II ~||2 • 
IMI2 

Moreover, C(p = \(p in the weak sense. 

Proof. As usual, in ( I3.7P we may restrict ourselves to the case H^lh = 1. Now, let 

A = {<p E W h2 (n) : y\\ 2 = 1, I ipdx = 0}. 

jSfep 1. First we show that J(ip,ip) is bounded from below on ^4. Repeating the 
estimates from the proof of Proposition 12.31 we obtain 

M / V/C(^)V^da; < -|| V^ 2 , + — 1| V/C|||^ L2 |M| 2 . 

Hence, for every <£> G A we have 

1 M 2 M 2 

j{<p,<p) > -||v^|| 2 - — liv/ciii^ii^n 2 > -— ||V/C||| 2 ^ 2 . 

Step 2. Let {y?})neN C A be a minimizing sequence that is 

A = lim J(<p n ,ip n ). 

n— too 

We show that <y9 n is bounded in W l,2 (Q). Since ip n is the minimizing sequence, there exists 
a constant C such that 

1 M 2 
C > J(<p n ,<p n ) > -||V^ n || 2 - — ||V/C||i 2 _ i2 , 

so we obtain 

\\V<p n \\ 2 2 <2C + M 2 \\V)C\\ 2 L ^ L2 . 
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Thus, using the Rellich compactness theorem we have a subsequence, again denoted by 
(p n , converging to (p strongly in L 2 (Q). Moreover, by the Banach-Alaoglu theorem, we 
obtain, again up to subsequence, also weak convergence of <p n towards to (p in W 1 ' 2 (Q). 

Notice, that (p G A. Indeed, by the weak convergence in W 1,2 {VL) we have that (p G 
W 1,2 (tt) and by the strong convergence in L 2 (Q) the limit function satisfy ||y5|| 2 = 1 and 
f n fidx = 0. 

Step 3. Now, we show that lim n ^. 0O J(<p n , <p n ) = J(<p, (p). 

First, notice that by the weak convergence of V</2 n in W 1,2 (Q) we have 

(3.8) liminf ||V^n||2 > ||V^|| 2 . 

n— >oo 

Next, by the strong convergence of tp n in L 2 (Q) and the fact that V/C : L 2 (Q) — > L 2 (Q) 
is linear and bounded it is easy to verify that 

V/C(ya„) — > V/C(<^) as n — )■ oo strongly in L 2 (Q). 

This property and again the weak convergence of ip n implies that 

V/C (<p n ) Vip n dx —7- / VK{(p)V(p dx as n — > oo 



which by estimate (13.81) together with previous step completes the proof of Step 3. 

Step 4- Finally, we show that the limit function (p satisfies the following eigenvalue 
problem L(p = \<p in the weak sense, namely 

J{(p,v) = A f fivdx for all v G W 1 ' 2 ^). 
Jn 

Let us denote 

J{(p + ev, (p + ev) 



f(t) 



f a (fi + ev) 2 dx 

for any v G W 1 ' 2 and eel This function is differentiable with respect to e near e = 
and has a minimum at 0. Hence the derivative vanishes at £ — 0, and we get 

n J (^) A<P,<P) L^ vdx T f~ \ \ f ~ A 



Hence the proof of Lemma 13.11 is finished. □ 
Now we are in the position to prove the Theorem 12.51 



Proof of Theorem \2.h\ As a standard practise, we show that under our assumptions, the 
linear operator C defined by the form f)3.6p has a negative eigenvalue A. Then, the function 



(Vwi(x)) 2 dx — M / / V x K(x, y)wx{y)Vw\{x) dy dx 
Jn Jq 

Ai / (wi(x)) 2 dx — MAi / / K(x,y)wi(y)wi(x) dy dx. 
Jq Jq Jq 
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ip(x,t) = e~ xt <p(x) with the eigenfunction if of £ corresponding to the eigenvalue A, is a 
solution of the linearized problem ( I2.8p -( 12TT0|) such that 

||^(-,t)|| 2 = e- A *||^|| 2 ^> oo. 

To do so, we use the definition of an eigenvalue of operator £ from Lemma 13.11 In view 
of ( 13. 7p . to prove that A < 0, it suffices to show that there exist ip G W 1,2 (Q) that 

J(cp,cp) < 0. 

Here, we choose y?(x) = Wi(x), where W\ is the eigenfunction of —A on Vt under the 
Neumann boundary condition satisfying j Q w\ dx = 1 and corresponding to the first non- 
zero eigenvalue Ai. Then, we obtain the following relation 

J(w 1 ,wi) = 



Now, since Ai > and j n (wi) 2 dx = 1, using assumption (12 . 1 3[) and choosing M > 1/A 
we complete the proof. 

□ 

4. Existence of solutions 

We construct local-in-time mild solutions of ( jl.ip - (jl.3p which are solutions of the fol- 
lowing integral equation 

(4.1) u(t) = e tA u - f Ve ( *" s)A (uVv) (s) ds 

Jo 

where e iA is the Neumann heat semigroup in fl. Moreover, we use the following estimates 
of {e iA } t > . 

Lemma 4.1. Let Xi > denote the first nonzero eigenvalue of —A in Q under Neumann 
boundary conditions. Then there exist constants C\, Ci independent of t, f which have 
the following properties. 

(i) If 1 < q < p < +oo then 

(4.2) ||e* A /|Up (n) <C(l + r 'U-f;)||/|| 

holds for all f G L 9 (fi). 

(ii) If 1 < q < p < +oo then 

(4.3) \\d x e tA f\\ LP{n) < cH(«-?)-3e-^||/|| M(n) 
is true for all f G L q (Q). 



STABILITY OF SOLUTIONS TO AGGREGATION EQUATION 13 

Proofs of above inequalities ( 14. 2 p and (14. 3 p are well-known and can be found e.g. in 

First, we construct global- in-time solutions in the case of mildly singular kernel. 

Proof of Theorem \2.1(K We split the proof into two parts. First we construct the local- 
in-time solution to problem (jl.ip -f fOj) and later on we show how to extend this solution 
on every time interval [0,T]. 

Step 1. Local-in-time solution. Here, we follow the reasoning from [121 Theorem 2.2]. 
We construct the local-in-time solution to the equation (14. ip . written as u{t) = e tA u + 
B(u,u)(t) with the bilinear form 

(4.4) B(u,v)(t) = - [ Ve MA («V/C(t)))(s) da, 

Jo 

in the space 3^t- Notice that e tA u G 3^r by ( I4.2p . To apply ideas from [T21 Theorem 2.2], 
one should prove the following estimates of the bilinear form (I4.4p . 
First, let us notice that by Minkowski's inequality we have that 

(4.5) 

\\V)C(v)\\ ql < || / \V x K(;y)\v(y)dy\\ q/ < [ \\V x K(;y)\\ q ,\v(y)\ dy < || V x Ji lU^MIx 
Jn Jn 

Now, for every y?, using (14.31) combined with relation (14. 5 p we obtain 

\\B(u,v)(t)\\x<C [ (t- s)- l/2 \\uV1C(v) (a) ||i da 
Jo 

<C [\t-s)- 1 ^u(s)\\ q \\V}C(v)(s)\\ q ,ds 
Jo 

< cuv^iu^ ['(t-sy^Wu^Uvis^ds. 

Jo 

Therefore, by the argument using in [12] we obtain 

(4.6) ||5(u, < CT^^-^HV^IIoo^llull^llull^ 

where |(1 - d(l - ±)) > 0. 

In a similar way, we prove the following L 9 -estimate 

(4.7) t^ l -^\\B{u,v){t)\\ q <Ct^ 1 -^ / {t-sY l/2 \\uVlC{v){s)\\ q ds 

Jo 

<C^ (1 -i) I {t-s)-V 2 \\u\\ q \\VK{v){s)\\„te 
Jo 

^eiv^iu^ 1 ^) /(f-^-^ii^iij^n.da 

Jo 



14 



RAFAL CELINSKI 



since 




V x K(-,y)v(y)dy\\ O0 < esssup\\V x K(x,-)\\ q '\\v 



Again, by the argument using in [12] we obtain 



(4.8) 



t^\\B(u,v)(t)\\ q < CT^ d{1 ^\\V x K 



■'\\ u \\yT\\ v \\yT- 



By inequalities (14.61) and (14. 8 p we obtain the following estimate of the bilinear form 



solution in y T by [121 Lemma 3.1]. 

Step 2. Global solution. Now, it suffices to follow a standard procedure which consists in 
applying repeatedly previous step to equation (II. ip supplemented with the initial datum 
u(x, kT) to obtain a unique solution on the interval [kT, (k + 1)T] for every fceM. Notice, 
that we can pass this procedure since the local existence time T depends only on ||u ||i 
and || V s -fr||oo,9' which implies that it does not change for all nonnegative no G L l (Q) with 
the same L 1 -norm (see Remark 12. ip . 



Now, we prove local-in-time existence of solutions in the case that K, is strongly singular. 

Proof of Theorem \2.11[ We assume now, that q' G [l,d]. Again notice that e tA uo G Xt 
since by (14.21) we have 



^^H^^CT^^-^IIV^ 




□ 



e* uq\\ Xt < C(||u ||i + ||tto|| 9 )- 



Next, for every u, v G 3^t, we get 



||B(u,t;)(f)||i<C / (t 



S )-^ 2 \\uVlC{v)(8)\\xd3 




< C\\V X K 




„ [ (t-s)- x ' 2 \\u{s)\\ q \\v(s)hte 



< CT 



/2 \\V X K 



Jo 

\oo,q'\\u\\x T \\v\\x T 



where C is a positive constant. 
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To deal with the L 9 -norm of B(u,v) we proceed similarly 

\\B(u,v)(t)\\ q <C /V-s)~ 1/2 y«V/C(i;)(s)||,da 
Jo 

< CWVvKWnj f (t-s)-^ 2 \\u(8)\\ q M8)\\ q ds 

Jo 

< cr 1/2 ||v a: -R'|| 00i(? '||M||A' T IK'IU T - 

Summing up these inequalities, we obtain the following estimate of the bilinear form 

||B(u,t;)||^ T < CT 1/2 ||V x .J^||oo, g '||M|U T ||w|U T - 

Hence, choosing T > such that 4CT 1 ' 2 ||V a; .fiT||oo,g'(ll' u o||i + ||uo||g) < 1; we obtain the 
solution in Xt by [121 Lemma 3.1]. □ 
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